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LONG-TERM  GOALS 

The  wave-driven  dynamics  of  the  coastal  ocean,  which  is  important  for  transport  processes,  mixing 
and  circulation,  is  strongly  affected  by  mud  deposits  on  the  continental  shelf  and  in  the  nearshore. 
However,  the  mechanics  of  wave-mud  interaction  and  the  impacts  on  coastal  wave  dynamics  is  not 
well  understood.  The  overall  objective  of  this  work  is  to  contribute  to  the  understanding  of  the  macro¬ 
scale  damping  effects  of  mud  on  waves  and  to  improve  modeling  and  prediction  of  wave  evolution 
along  muddy  coastlines. 

OBJECTIVES 

The  specific  objectives  of  the  proposed  effort  are  to:  1 )  establish  a  comprehensive  dataset  from  two 
new  datasets  acquired  in  2008  (ONR-funded  wave-mud  MURI)  of  waves  propagating  across  more 
than  25  km  over  a  muddy  seafloor,  and  2)  determine  the  characteristics  of  mud-induced  dissipation 
of  wave  energy  from  intermediate-  to  shallow-water  depths. 

APPROACH 

To  improve  our  understanding  of  mud-induced  damping  rates  on  surface  waves,  and  its  variation 
across  the  shelf,  we  combine  observations  made  during  the  2008  MURI  field  experiment 
(MUDEX08)  by  the  NPS/SIO  and  WHOI  teams,  to  establish  a  comprehensive  data  set  of  wave 
evolution.  The  high  spatial  coverage  allows  a  detailed  analysis  of  mud-induced  wave  damping 
across  the  shelf  to  the  shoreline,  including  the  transition  to  shallow  water.  The  proposed  work 
includes:  1)  a  detailed  analysis  of  the  combined  data  set  to  determine  damping  characteristics 
across  the  shelf,  and  2)  implementation  of  a  conventional  third-generation  wind-wave  model  to 
assess  the  effects  of  mud  on  wave  propagation  and  wind-wave  generation,  and  3)  investigate  the 
role  played  by  nonlinear  interactions  or  other  coastal  processes  in  the  damping  of  short  waves. 
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Figure  1.  (a)  The  study  area  is  located  west  of  the  Atchafalaya-Vermillion  Bay  system  in  the  Gulf  of  Mexico, 
(b)  Bathymetry  (black  curves  are  isobaths ,  units  m)  and  sensor  locations.  Blue  dots  are  the  inner  shelf  stations 
where  dw-stations  are  Datawell Directional  Waverider  buoys,  the pv  are  Nortek  Vector  ADV-pressure  sensors, 
and  pa  are  pressure  sensors.  The  nearshore  array  of  SONTEK  Triton  ADV-pressure  sensors  is  indicated  with 
red  squares.  The  green  triangle  shows  the  location  of  the  meteorological  buoy. 
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WORK  COMPLETED 

Combining  the  2008  MURI  data  sets 

The  datasets  collected  by  the  NPS/SIO  and  WHOI  teams  during  the  2008  MURI  field  experiment  on 
the  Louisiana  shelf  (see  also  Trainor,  2009;  Trainor  et  al.  2008;  Garcia-Garcia  et  al,  2012),  were 
combined  into  a  comprehensive  dataset.  The  combined  data  set  includes  a  three-transect  array  on 
the  inner  shelf  (Figure  1,  inset  1)  and  a  high-resolution  nearshore  array  that  smoothly  connects  to 
the  shallow-end  of  the  western  transect  of  the  shelf  array  (Figure  1 ,  inset  2).  The  new  data  set 
includes  observations  with  a  high  spatial  resolution  (see  figure  2)  across  many  wavelengths 
(approximately  25  km)  and  will  provide  a  much-needed  database  for  the  study  of  the  evolution  of 
waves  over  mud.  The  dataset  consists  of  observations  of  a  wide  range  of  wave  and  wind  conditions 
(figure  3),  including  fetch-limited  wave  growth  (wind  from  northerly  directions),  swell  propagation 
(southerly  waves,  weak  winds),  and  mixed  sea-swell  events. 


Figure  2  Boundary  conditions  for  the  southern  boundary  (red  line)  were  taken  from  observations  (dw12  and 
pv16).  1 D  non-stationary  runs,  initiated  by  dw12  and  pv16  for  the  western  and  eastern  boundary  (blue  lines), 
respectively,  were  used  as  side-boundaries. 

Implementation  third-generation  wave  model 

To  study  the  effects  of  mud  on  the  seafloor  on  the  nearshore  wave  energy  balance  in  more  detail, 
we  have  implemented  a  conventional  third-generation  wind-wave  model  (SWAN).  Wave  boundary 
conditions  are  taken  from  the  most  offshore  buoy  (DW12),  wind  forcing  is  obtained  from 
meteorological  observations  made  available  by  the  WHOI  team  (Drs  Trowbridge  and  Fredericks), 
and  bathymetry  information  was  taken  from  the  NOS  coastal  relief  model  augmented  with  nearshore 
observations  by  the  WHOI  team  (Elgar  2009,  personal  communication)  during  the  experiment. 

To  account  for  the  down-wind  variability  of  the  atmospheric  boundary  layer  due  to  the  decrease  in 
roughness  length  over  water,  wind  speeds  during  offshore  wind  events  (defined  as  wind  events  with 


mean  wind  directions  <  ±  —  from  exactly  offshore)  are  modified  by  a  spatially  varying  scaling  factor 

(Taylor  &  Lee,  1984).  The  model  was  run  in  third-generation  mode  (GEN  3)  with  saturation-based 
whitecapping  (Van  der  Westhuysen  et  al.,  2007)  combined  with  the  Yan  wind  input  term  (Yan,  1987). 
All  available  source  terms  are  included  in  the  computations  except  the  triad  interactions. 

The  present  model  implementation  was  used  as  a  reference  for  an  equivalent  sandy  shelf  to  identify 
the  principal  effects  of  the  mud  on  the  observed  wave  dynamics.  Thereto  we  ran  the  model  with  a 
standard  bottom  friction  term  (Hasselmann  et  al.,  1973)  with  the  (fixed)  friction  coefficient  set  to 
0.038  m2s'3  (Van  Vledder  et  al.  201 1 ). 

Analysis  of  cross-shelf  wave  dissipation 

To  analyze  the  effects  of  mud  on  the  wave  energy  balance,  we  defined  a  cross-shore  direction  at 
about  10°  clockwise  from  North  (see  figure  2),  and  assumed  straight  and  parallel  depth  contours  in 
the  alongshore  direction  so  that,  for  a  stationary  wave  field,  the  wave  energy  balance  reduces  to 

—  6)  =  e>gen(w,  #)  +  ^disO^S  6)  +  <S>ni(c<;,  6)  (1) 

where  T{lo,  9)  is  the  cross-shore  wave  energy  flux  and  the  forcing  terms  on  the  right  account  for 
(from  left  to  right):  wind  generation  (input),  dissipation  (whitecapping,  wave-bottom  interaction  etc.), 
and  nonlinear  interactions  (redistribution).  Through  comparison  of  the  observed  and  modeled  cross¬ 
shore  energy  flux  gradient  (1),  and  normalized  growth  rate,  we  assessed  the  macro-scale  effects  of 
wave-mud  interaction  on  evolution  of  waves  across  a  muddy  continental  shelf. 

We  have  further  considered  effects  of  the  nearby  Atchafalaya  Bay  system  and  the  sediment  plume 
dynamics  in  the  coastal  zone  on  the  observed  wave  damping  and  development. 
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RESULTS 


The  model-data  comparison  shows  that  wave  heights  are  generally  predicted  well  with  the  standard 
Jonswap  bottom  friction  term,  apart  from  episodic  events  when  dissipation  levels  in  the  observations 
appear  to  be  enhanced  (see  figure  3).  This  suggests  that  a  simple  (fixed)  friction  coefficient  is 
inadequate  to  capture  the  more  dynamic  (and  time  varying)  rheology  of  the  muddy  seafloor,  which  is 
consistent  with  previous  findings  by  other  researchers  (e.g.  Sheremet  et  al.,  2005).  Further,  during 
swell-type  events  (longer-period  waves)  the  muddy  seafloor  appears  to  affect  wave-bottom 
interaction,  and  results  in  time-varying  enhancement  of  losses  of  wave  energy.  During  times  that  the 
wind  is  blowing  in  offshore  direction  (fetch-limited  and  slanted  fetch  wave  growth)  the  effect  of  mud 
appears  to  be  a  suppression  of  wave  growth,  either  through  the  damping  of  short  waves  or  through 
decreasing  the  effectiveness  of  the  momentum  transfer  from  the  wind  to  the  water. 


Figure  3.  Significant  wave  height  versus  time  across  (a)  the  western  transect  in  1 1 .3-  (pv2,  red  curve),  8.3- 
(pv4,  black),  3.9-  (n15,  blue),  and  1.7-m  (n4,  green)  water  depths,  (b)  the  central  transect  in  10.9-  (pv7,  red) 
and  8.3-m  (pv9,  black)  water  depths,  and  (c)  the  eastern  transect  in  10.9-  (dw12,  red)  and  5.5-m  (pv16,  black) 
water  depths.  The  dashed  black  lines  at  1m  height  are  for  reference.  Data  gaps  for  shelf  stations  between 
March  2  and  March  5  are  due  to  instrument  maintenance.  Shaded  areas  refer  to  events  discussed  in  the  text 
and  indicate  periods  with  large  spatial  wave  variability  and  strong  (>  9  m/s)  winds  from  the  northeast  (yellow,  0 
<  dir  <  80  deg)  and  northwest  (gray,  300  <  dir  <  360  deg),  as  well  as  a  period  of  large  swell  (blue). 

The  model  hindcasts  are  in  fairly  good  agreement  with  the  observed  wave  height  variability  during 
the  experiment  for  all  stations.  However,  for  fetch-limited  conditions  (northerly  winds),  the  model 
tends  to  overestimate  wave  heights  (figure  4b),  whereas  during  onshore  wave  propagation 
(southerly  winds),  the  agreement  is  considerably  better  (figure  4a). 
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Hs  observations  [m]  Hs  observations  [m] 


Figure  4.  Modeled  versus  observed  significant  wave  height  for  all  stations  on  the  western  and  central 
transects  for  (a)  onshore  winds  (130°-250°  true  north)  and  (b)  offshore  winds  (310°  -70°).  (Cross-shore  is 
rotated  10°  clockwise  from  true  north).  The  dashed  black  lines  indicate  perfect  agreement.  The  slope  of  the 
best-fit  line  is  0.94  for  (a)  and  1.18  for  (b),  r2  =  0.94  for  (a)  and  r2  =0.71  for  (b),  and  the  root  mean  square  error 
is  0.11  m  for  (a)  and  0.19  m  for  (b). 

To  investigate  a  possible  correlation  of  coastal  sediment  plume  dynamics  with  the  observed 
differences  in  modeled  and  observed  wave  growth  at  higher  frequencies,  satellite  images  from 
MODIS  Terra  250  were  analyzed.  The  images  were  obtained  from  the  NASA  EOS  Data  Gateway, 
processed  with  HDFLook,  and  converted  from  percentage  reflectance  to  an  estimate  of  total 
suspended  matter  in  the  surface  layer  (Miller  and  McKee,  2004). 

From  the  available  satellite  data  (cloud  cover  limits  visibility),  it  is  seen  (see  e.g.  figure  5)  that  the 
westward  extent  of  the  plume  is  pushed  farther  west  during  easterly  wind  conditions,  resulting  in 
higher  surface  sediment  concentrations  at  the  experiment  site.  When  wind-driven  currents  carry 
higher  sediment  concentrations  westward,  wave  growth  owing  to  winds  appears  to  be  suppressed. 

In  contrast,  sediment  concentrations  are  relatively  low  at  the  experiment  site  (e.g.  March  8  (not 
shown),  and  February  28,  figure  5c)  when  modeled  and  observed  wave  evolution  are  in  better 
agreement,  including  in  the  higher-frequency  energy  balance.  Moreover,  satellite  imagery  (not 
shown)  suggests  that  the  location  of  the  western  sensor  transect  often  coincides  with  the  maximum 
western  extent  of  the  Atchafalaya  sediment  plume  on  the  shelf  (the  plume  extends  farther  westward 
closer  to  shore).  In  addition,  the  western  transect  also  corresponds  to  the  western  most  extent  of 
Atchafalaya  sediment  deposition  on  the  shelf  (Draut  et  al.  2005),  possibly  explaining  the  rapid 
“cleaning”  of  the  surface  waters  in  the  study  area  during  westerly  winds,  and  therefore  the  sensitivity 
to  changes  in  the  wind  direction  in  the  model-data  comparisons. 
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Figure  5.  Estimate  of  total  suspended-matter  concentrations  (from  Modis  250  imagery)  on  (a)  March  15,  (b) 

March  16,  and  (c)  February  28.  All  images  were  recorded  around  17:00  o’clock  GMT.  High  concentrations  are 
red,  low  concentrations  are  blue  (color  scale  on  right).  Grey  shading  is  land.  Instrument  locations  are  shown 
with  black  circles  (shelf  stations)  and  red  squares  (shoreward  stations). 

The  frequency  distribution  of  the  dissipation  (see  figure  6)  between  observations  and  model  shows 
good  agreement,  although  the  model  underestimates  dissipation  toward  higher  frequencies  (>  0.2 
Hz)  at  the  inner  shelf  stations.  It  can  be  seen  (figure  6)  that  dissipation  takes  place  in  the  energetic 
ranges  of  the  spectrum.  The  distribution  of  dissipation  is  mostly  consistent  with  a  direct  interaction 
mechanism  where  damping  is  induced  through  interaction  of  wave-induced  fluid  motions  with  the 
seafloor,  such  as  is  the  case  for  bottom  friction  (Hasselmann  et  al,  1973),  or  a  direct-interaction  two- 
layer  model  (e.g.Dalrymple  &  Liu,  1978;  Ng,  2000). 
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Figure  6.  Contours  (color  scales  on  the  right)  of  (a  and  d)  observed  energy  flux  gradients,  (b  and  e)  modeled 
(SWAN)  energy  flux  gradients,  and  (c  and  f)  observed  energy  density  as  a  function  of  frequency  and  time,  a-c 
are  between  1 1 .3  and  8.9  m  depth  and  d-f  are  between  3.9  and  2.5  m  depth.  Positive  flux  gradients  indicate 
generation,  negative  flux  gradients  indicate  dissipation.  Note  the  different  scales  of  (a)  and  (b)  versus  (d)  and 
(e). 
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IMPACT/APPLICATIONS 

The  availability  and  analysis  of  a  comprehensive  community  dataset  of  wave  propagation  across  a 
muddy  shelf  will  be  an  important  contribution  to  the  study  of  wave-mud  interactions  and  the 
validation  and  calibration  of  new  theories  and  modeling  approaches. 

The  fact  that  the  time-varying  properties  of  mud  are  important  even  for  the  macro-scale  effects  of 
wave  propagation  in  coastal  areas,  and  cannot  be  adequately  represented  by  a  fixed  rheology,  is 
important  for  operational  wave  prediction  in  coastal  areas. 
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